The Behavior of Electronic Interferometers in the Non-Linear Regime 
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We investigate theoretically the behavior of the current oscillations in an electronic Mach-Zehnder 
interferometer (MZI) as a function of its source bias. Recently, The MZI interference visibility 
showed an unexplained lobe pattern behavior with a peculiar phase rigidity. Moreover, the effect 
did not depend on the MZI paths difference. We argue that these effects may be a new many-body 
manifestation of particle-wave duality of quantum mechanics. When biasing the interferometer 
sources beyond the linear response regime, quantum shot-noise (a particle phenomena) must affect 
the interference pattern of the electrons that creates it, as a result from a simple invariance argument. 
An approximate solution of the interacting Hamiltonian indeed shows that the interference visibility 
has a lobe pattern with applied bias with a period proportional to the average path length and 
independent of the paths difference, together with a phase rigidity. 
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In the last two decades, electron interferometers be- 
came a primary tool in mesoscopic physics, for in- 
vestigating quantum coherence of transport in semi- 
conductors, and to measure and control novel quantum 
effects [H, H, S S 0|- Their behavior is understood 
so far only in the linear response regime, where a small 
bias is put on their sources, and the non-interacting pic- 
ture is valid. In the non-linear regime several electrons 
are present inside the interferometer at a given time 
and may form, in contrast to the optical interferome- 
ters, non-trivial many-body correlations due to Coulomb 
interaction. Indeed, new experiments in the non-linear 
regime demonstrated that the Landauer-Buttiker formal- 
ism Q seems to break down, and the interference pattern 
showed new peculiar behavior. Recently, an unexpected 
interference behavior of a Mach-Zehnder interferometer 
(MZI) was reported in which the visibility (propor- 

tional to the observed amplitude of the Aharonov-Bohm 
(AB) oscillations of the drains current) evolved in a lobe 
pattern with increasing the source bias, with zero visibil- 
ity between the lobes, and a phase independence on the 
bias inside each lobe ("phase rigidity"). Similar visibil- 
ity lobes were observed earlier in a two-terminal closed 
interferometers [2, In the MZI case, the lobe pattern 
could not be explained using any non-interacting picture, 
mainly because the lobes were very robust against in- 
duced asymmetry between the two path, and did not 
depend at all on their length difference. Recent expla- 
s proposed : 

that were geometry dependent, somewhat in contrast to 
the robustness of the experimental observation. More- 
over, these propositions do not explain the total phase 
rigidity seen in the experiment Q. 

Here we show that these lobes and phase rigidity are 
indeed very reasonably a non-perturbative result of the 
Coulomb interaction. We argue that in the non-linear 
regime, this interaction causes the extra electrons in- 



nations proposed for this effect using Bosonization tech- 
niques [H , 13, lH imposed interactions between electrons 



side a two path interferometer to form a correlated state 
which is rooted in their quantum behavior; as electrons 
behaves both as particles and as waves, the phase of an 
interfering electron gets discrete "quantum kicks" accord- 
ing to the possible occupations of the additional electron 
states in the two arms. In some distinct source voltages 
those phase fluctuations cause the vanishing of the in- 
terference, which reappear at lower or higher voltages, 
leading to a lobe pattern in the visibility. The basis for 
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FIG. 1: Schematic of the MZI geometry, with generally asym- 
metric two paths. 

our derivation is a general principle for mesoscopic de- 
vices, which we term Buttiker-Levinson argument (BLA): 
changing the chemical potential of all the sources by the 
same amount, should not change the transmission from 
any source to any drain. It is generally relevant for any 
energy-dependant mesoscopic devices [141], and in partic- 
ular for the MZI, being two-path interferometer formed 
by one dimensional edge channels in the quantum Hall 
effect regime, since the interaction between the electrons 
inside the paths is unscreened by the gates nearby. Sur- 
prisingly, when applied to electron interferometers in the 
non-linear regime, this simple argument is violated by 
non-interacting models Q. 

A seemingly natural solution to this inconsistency is to 
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correct the non-interacting models by taking into account 
a mean field approximation (MFA)@,[l3,[y| in which the 
bias induces a static mean charging in each of the interfer- 
ometer arms which causes a phase shift due to Coulomb 
interaction. Experimentally however, this approximation 
failed for the MZI as the lobes were independent on the 
two arm length difference [1] . We now further argue that 
this MFA correction is never valid in the MZI, since when- 
ever it is significant, quantum fluctuations are strong and 
cannot be neglected. As the phase shift must restore the 
BLA validity, the phase shift correction in the MZI must 
be given by e i 2n (Ni-N 2 ) jjj^ w here N a is the excess num- 
ber of electrons in path a due to the source bias. Due to 
the applied bias this number of electrons and the phase 
in each arm fluctuate binomially (quantum shot-noise) 
with a large variance and, most importantly, with strong 
negative correlation between the two arms. Therefore 
the induced phase shift cannot be approximated simply 
by the average charging since ^e tS ^ N ') ^ ^(^W) ^ as 
was suggested in [lflj. Rather, in this situation of few 
electrons inside the MZI generally [N, S<fi] 7^ 0, and non- 
linear interactions may cause new many-body coherent 
states analogous to squeezed spin states 0. Our sit- 
uation is shown below to be closely related to another 
non-linear system, a biased finite ID detector channel 



coupled a MZI arm [17|, |1 



Our starting point is the non-interacting Hamiltonian, 
written for spinless or spin polarized electrons in the MZI 
drawn in Fig. [T] It is expressed by the creation and anni- 



hilation operators c 



for the incoming extended 



single-particle energy states from the two sources £1,2, 



(S Q )t (S a ) 



(1) 



a=l,2 k 



In this model a linear relation can be established between 
the incoming and outgoing fc-operators from the sources 
and to the drains, respectively 
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where the s-matrix for the particular MZI geometry can 
be calculated from the product S ( MZI *> = S (Q PC2 ) • ■ 
s (QPCi) w fth the three s-matrices defined as (a = 1, 2) 
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where r a and t a the reflection and transmission amplitude 
of QPC(a), and L\^. are the two MZI path lengths. We 
can define the annihilation operator c£ of an electron 
moving in arm a of the MZI as obeying the relations 
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Considering the local nature of the interactions and of 
the device, we choose to express the current in the drains 
using the local operators after QPC2, ipm a <\(x,t) = 



(t). The resulting chemical potential at any 



drain D a , /iD a , is proportional to the average current to 
the drain, fip, a = -(l£> a ), and since our problem is sta- 
tionary in time it can be written (linearizing e(k)) using 
this operator ip/p>\{x,t) at t = 



fJ-D Q 
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where v g is the group velocity at the Fermi energy. In 
the non-interacting model Eq. [5] can be written as 

M-Dc, = Mo,a + ^(AV,a eI¥, °)' with Mo, a a real parame- 
ter and n v , a a complex pre-factor of the phase (ipo = 
kpAL) dependant term. The second term is proportional 

to J^cos^o + ^(e-E F ))de = [viv.dn^)] " 
cos^o + ^t(m — Bp)), with Ep the Fermi energy at zero 
bias. It violates BLA through the explicit dependancc 
on the average source bias /2 — Ep. 

Since, as stated above, a MFA correction is inadequate 
due to the quantum fluctuations, we are forced to return 
to the Hamiltonian, and introduce additional non-linear 
Coulomb terms (whose range we assume here to be longer 
then the MZI arms) to take care of the BLA 
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where N a (t) = J Q a p a (x,t)dx is the number operator 
which counts the electrons which are added to the MZI 
path a region at time t, p a (x,t) — ipl l (x,i)tl> a (x,i), and 
Ca = - rTr 2L is the electric capacitance of path a, induced 

flVg 

by the dispersion in Hq and Pauli principle. Hence, the 
added terms in the Hamiltonian do not have additional 
free parameters. 

The role of the additional term is to locally raise the 
bottom of the conduction band (the electro-static poten- 
tial) in path a according to its overall local charging. 
The effect is exactly such that when a full beam between 
energies Ep and Ep + A/i enters the path, the last fully 
occupied state remains with same momentum kp, with a 
new energy p, = Ep + A/i, which restores BLA validity. 
We should note that in addition to this dynamic effect, 
the new terms also add electron-hole excitations above 
the Fermi sea, which we ignore [l9j |. 

We shall solve now approximately the equation of mo- 
tion (EOM) for the interacting case. Keeping the defini- 
tion in Eq. [5] and expecting the same functional depen- 
dence of jiB a on ipo as in the non-interacting case, we 
want to show that (i v , a develops a lobe structure behav- 
ior and phase rigidity with an energy scale independent 
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of AL. The EOM for ip a (x,t) reads 



operator 



di/i a (x,t) 1 

= —[1pa{x),H] = 
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dt ih 
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Where <d(x) is the Heaviside function. We use the linear 
relations in Eq. 3] to define the operators ip a (x,t) at 
x < and at x > L a 

V'CDjOM) = X! S «| P V^OM): ^>io (8) 
(8=1,2 

V>(co(M) = S % PC ^(SpjOM)) x<0 

0=1,2 

where ip(s a )(x,t) = '52 k e. lkx c'^''\t). Having the non- 
interacting solution ip a (x,t) = ipa\x — v g t) at x < as 
boundary condition, Eq. [7] has a unique formal solution, 
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Mx,t) = { T[U a (t)]V& 0) (x -Vg^T- 1 \U a {t)\ 0<x<L o 
if) a (L a ,t - (x - L a )/v g ) L a < x 

where T is time-ordering operator, and we define 



U a (t) = e 



(9) 
(10) 



Eq. [5] is difficult to evaluate analytically, due to the fact 
that U a (t) depends on ip a (x, t). However, one can obtain 
some general properties of the highly correlated quantum 
state at the drains, using the following derivation. First 
we approximate U a (t) by replacing p a (x, t) in N a (t) with 
the non- interacting density p® a (x — v g t). One should note 
that this is a very crude approximation as electrons really 
repel each other and change each others group velocity 
considerately; hance the true solution might be different 
in its fine details. However, note that unlike the MFA, 
this approximation maintains the non-Gaussian behavior 
of the shot-noise, as well as the negative correlation of the 
shot-noise between the two paths. Next, due to the fact 
that N°(t) are Bosonic, [N°(t), N°(t')] is a C-number, 
so the T-ordering can be eliminated by repeatedly using 
Baker-Hausdorff formula, such that the phases originat- 
ing from the T and the T _1 ordering cancel each other. 
Then, we take t = and rearrange the operator in the 
exponent U a = U a (t = 0) = e~ t ^ a as a weighted inte- 
gral over the density operator <p a = J_ w a (x)p° a (x)dx. 
Straightforward calculation shows w a (x) to be a triangle 
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An additional approximation is to restrict the energy ex- 
citations in the operators p° a (x) and therefore also in the 
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where w a ^-k' is the Fourier transform of w a (x), and O 
denote an operator restricted to the energies [Ep,Ep + 
Ap] . This approximation is very intuitive as it takes into 
account only dephasing due to excess shot-noise, disre- 
garding virtual transitions to energies above Ep + Ap. 
It was proven very useful in explaining the unique ex- 
perimental results of dephasing by shot-noise 17 1. With 
these approximations, the pre-factor of the phase depen- 
dant part of Eq. [5] reads 

AV = hvg(y \uiu 2 ^ (£i) $1 (L 2 ) UlU 3 \*o) (13) 



where |* ) = (jltt^* c< a ») fo-CWx,, - E F )). 

Eq. [13] can be evaluated by noting that 4>i and 4>2 com- 
mute, so one can introduce a basis of a complete set of 
eigenstates of the two operators {|y>i, ¥2)}- Note that the 
single-particle matrix w a ,k,k' , restricted to the voltage 
window [Ep, Ep+Ap], has a discrete set of non- vanishing 
eigenvalues, corresponding to localized electron-states 
near the finite influence region of w(x) (Eq. [TT|) . The 
eigenstate \ip a ) of the operator cj> a is a specific choice of 
occupations of those single-particle states, the value ip a 
being the sum of the eigenvalues of the occupied states. 
With two insertions of such a full set, Eq. [TBI now reads 

AV = hv g (*o|fPi,^>(vi|^? + (ii,0)|v»i)(14) 

Vl;V2;V' 1 ;V' 2 

■fal® (£2,0) |^><¥>i,^|*o) • e^-^'^ 

This expression differs from the non-interacting one only 
by the exponent at the end of the r.h.s. of the equation. 
Eq. [14] can be evaluated numerically, through diagonal- 
ization of the matrices w a ,k,k' ■ Fig. [2]shows the result for 
the visibility, as a function of (N) = (N 1 +N 2 ) = j^Ap, 

at = TT = °- 2 > and l*i I 2 = 1*2 1 2 = 0.5. The pre- 

diction of the non-interacting model is also plotted for 
comparison. One can clearly see that while in the non- 
interacting model the visibility and phase are slowly vary- 
ing (because the scale on which they change is propor- 
tional to AL -1 ), in the interacting model a lobe pattern 
appears in the visibility, with a stick-slip behavior of the 
phase. The lobes evolution is much faster, apparently 
proportional to (N)^ 1 and hence to L , in agreement 
with the experimental results. Fig. [3] shows the evolu- 
tion of the eigenvalues of the single-particle matrices W\ 
and W2 as a function of ( N) . At small bias there are two 
non-zero eigenvalues, one for each of the phase opera- 
tors. When put in the exponents of Eq. Q~U they exactly 
cancel the kinetic phase induced from Hq, resulting in 
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FIG. 2: The amplitude and phase shift of the interference 
oscillations in the chemical potential of the MZI drain, as 
a function of (TV) = LA/i/hvg in the non-interacting model 
(red curves) and in presence of the interaction (Eq. (|14|l , black 
curves), for (AN) /(TV) = AL/L = 0.2. The lines between the 
calculated points are a guide to the eye. Note that we plotted 
here the total current visibility, which is more suitable for 
analyzing the lobes then the differential visibility [IJ. 
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FIG. 3: The pairs of the single-particle eigenvalues of the 
phase operators wi and W2 of the two paths (blue and green 
lines), as a function of (TV), for (AN)/{N) = 0.2. Black lines 
mark the zero visibility in Fig. ((2). Note that at every such 
zero, there is a pair of eigenvalues whose sum is 7t/2 



phase rigidity. The first zero of the visibility occurs at 
(TV) = 1. As seen from Fig. [3l at this value the sum 
of the second pair of single-particle eigenvalues of the 
two arms is exactly X\ + A, = tt/2. This can be un- 
derstood, by assuming that the operators i//^ (Li, 0) and 
■02 (L2, 0) in Eq. [I4]act effectively only on the first single- 
particle state (the numeric calculation indeed shows that 
the contribution of the other electron transitions in Eq. 
[T4l is negligible) . However the phase of these two oper- 



which 



ators is modulated by tp^ — ip^ 
contains also the occupation of the second eigenstates. 
As the second electron can occupy the lower arm or the 
upper arm, a zero visibility happens when the difference 
in the added phase between these two options become 



7r, so they coherently cancel each other. Quite generally, 
the visibility goes to zero whenever a sum of a pair of a 
single-particle eigenvalues of the two paths reaches ir/2. 

Compared with the experiment, the theory still has 
shortcomings. The first is that the visibility lobes do 
not have the overall decaying envelope. Moreover, dif- 
ferentiating the curve in Fig. ((2]) near (TV) = 1, leads 
to a phase-dependant differential response d^iD a ,ip/dfJ.s 1 
which is twice that of the linear response; a differential 
visibility of 200%. These might be a nonphysical result 
of our approximations, which will be absent in the exact 
solution of Eq. 

In conclusion, in this paper we investigated the effect of 
quantum shot-noise inside electronic interferometers on 
the interference visibility in the non-linear regime. It was 
well established that an interaction Hamiltonian such as 
in Eq. [6] must always replace the non-interacting one, 
which quite generally leads to lobe pattern of the visibil- 
ity with increasing source bias. It would be desirable to 
apply the theory developed here to other mesoscopic de- 
vices such as the two-terminal AB interferometers [1, , 
the MZI working in the fractional quantum Hall regime, 
or the two-particle interferometer p, @] . 
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